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ABSTRACT 

A simple system was designed to collect wide-angle reflection 
records in order to investigate the interval sound speeds of the sedi- 
ment layers in the ocean. The system consisted of a frequency 
modulated receiver, a cut -to -channel Yagi antenna, and a sonobuoy, 
used in conjunction with a precision sonic profile recorder and a 
triggered sound source. A computer routine for reducing the data was 
obtained and modified for compatibility with the system and use on the 
IBM 360/67 computer. The system was designed to be both inexpensive 
and simple to use without any loss of accuracy. 
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I. INTRODUCTION 



The objective of this study was to develop a simple yet in- 
expensive system for collecting wide-angle reflection data . This 
system would be used in conjunction with normal-reflection profiling 
equipment normally carried on board oceanographic survey vessels. 

The study involved testing an inexpensive radio receiver to be used to 
receive a signal transmitted by a sonobuoy with enough clarity that 
the signal could be recorded on a precision seismic graphic recorder. 
The study also involved testing and launching five sonobuoys and in- 
vestigating modifications which can be made to the sonobuoys to 
enhance the success of the system. 

A further objective of the study was to provide a method for re- 
ducing the data collected by the system and obtaining the mean sound 
speed existing between successive reflectors (called the interval 
sound speed of the layer defined by an upper and lower reflector) and 
the thicknesses of the corresponding sub-bottom layers. This required 
modifying an existing program to read the data and compute the results 
on the IBM 360/67 digital computer using FORTRAN IV language. 

This paper presents a short discussion of the background and 
purpose for developing a system for collecting wide-angle reflection 
data and obtaining a program to reduce the data and determine the in- 
terval sound speeds. A theory of wide-angle reflections is offered, 
together with a mathematical model used to develop the equations and 
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suggest a method of solution. The computer program is discussed, 
including the method of digitizing data from a wide-angle reflection 
record. The radio receiver, antenna, and radio sonobuoy used to 
develop the system for use in conjunction with normal-reflection 
equipment is discussed. In the last section, conclusions and sug- 
gestions for further work on the subject are offered. 
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II . BACKGROUND 



The geological structure and composition of the sub-bottom 
layers in and around Monterey Bay are presently known to only a 
limited extent. Sufficient depth recordings have been made to enable 
drawing a detailed relief map of the sea floor; and shallow cores and 
grab samples have been analyzed for determination of the composition 
of the sediments in the first few meters . But knowledge of the deeper 
layers has been restricted to only subjective analysis involving the 
extrapolation of the immediate terrestrial structure from wells and out- 
crops and correlating the data with normal-reflection seismic profiling 
data. Unfortunately, normal-reflection profiles measure only the 
travel times and relative intensity of the vertical reflections. Without 
knowledge of the sound speed-depth relationship, only a guess can be 
offered as to the actual sub-bottom structure. A knowledge of the 
sound speed-depth variation is essential to a quantitative evaluation 
of the sediment layers and to a geological interpretation relating sound 
speed to type of material, age, and consolidation (Le Pichon, Ewing, 
and Houtz , 1968) . 

The system developed for collecting wide-angle reflection data 
along with the computer program to reduce the data can provide in- 
formation on the sound speed-depth relationships and the layer thick- 
nesses . 
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III. THE WIDE-ANGLE REFLECTION METHOD 



The characteristic feature of the wide-angle reflection method is 
the measurement of the travel times of longitudinal waves which have 
been reflected at boundaries separating media of different acoustic 
impedances (the product of sound speed and density) . From measure- 
ments of reflection times it is usually possible to determine the depths 
and dips of the reflection horizons and the speed of the seismic wave 
(Jakosky, 1949). 

A. THE THEORY 

The system was developed based on the theory that a sound pro- 
duced at source S (Tig. 1) will reach the hydrophone of the sonobuoy B 
by three different ray paths as shown; J_ the direct ray which travels 
just below the surface of the water, 2_ the sound reflected from the 
water-sediment interface, and 3_ the sound waves reflected from the 
various sediment layers that may exist below the ocean floor. If the 
distance separating the sound source and the hydrophones is steadily 
increased while the source emits sharp sound pulses at regular intervals, 
the time difference between the arrivals of the direct and reflected 
sound waves will change depending on the sound speed -depth variation 
within the sediments. A record of time differences is used to compute 
the thicknesses and interval sound speeds of the various layers and 
obtain the sound speed-depth relationship. 
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Figure 1. Wide-angle reflections 
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B. 



THE MATHEMATICAL MODEL 



1 . Forming the T/X Plot 

For the purpose of forming an elementary concept of the 
reflection method and to illustrate the formation of the wide-angle re- 
flection records, it is useful to first consider a simple case of a 
single layer of homogeneous isotropic material of uniform thickness H. 

It follows from this case that the ray paths will follow straight lines. 

A further assumption may be made that the lower surface reflects the 
sound wave in accordance with the physical laws of optics so that the 
angle of reflection is equal to the angle of incidence. Also it will be 
assumed that the path taken by sound passing from one media to another 
having a different sound speed is refracted in accordance with Snell's 
Law . 



Figure 2 illustrates a single layer model with a water 

medium of depth H and average vertical sound speed V. The theoretical 

travel time on the record for a vertically-propagating sound wave when 

the sound source and sonobuoy are at the same point S in the ocean is 

referred to as T . Thus T defines the origin of the record at which 
o o 

point the direct travel time is zero. The source sound wave will 
propagate downward through the water to point A and then be reflected 
back to the sonobuoy in time 



T = 
o 



2H 

V 



( 1 ) 
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Figure 2 . Ray paths and T/X. curves for rays 
reflected from a horizontal surface 
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If an anology is made to the laws of geometric optics, the 



ray paths may be conveniently handled by utilizing the concept of 
images . The travel time to the theoretical image point O' is equal to 
the time for the reflected ray to return to the sonobuoy. Therefore, 
when a mathematical model is constructed, the reflected times can be 
considered as the time required for the sound to travel from any shot 

point S. to the image point O' at a sound speed of V. 

1 

At shot point S^, a sound wave travels to point where 
it is reflected upward to the sonobuoy. This wave is referred to as 
the reflected ray from the first reflector. The direct ray travels hori- 
zontally from the source to the sonobuoy at S. Again using the concept 
of images, the reflected ray path S^B^S is equal to the assumed path 
S^B^O' (path length R^) and the travel time becomes 



t i = 



R, 



V 



( 2 ) 



where T^ is the travel time of a reflected ray. From the geometry of 
the triangle SB^A 




2H 

cosj2fj 



(3) 
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Equation (2) then becomes 



T 



2H 



(4) 



1 VcosjE^ 



where 0 is the angle of incidence measured from the vertical. 

The direct ray arrives at the sonobuoy at time 



X 



1 



D 



1 



V 



(5) 



where is the travel time of a direct ray, X^ is the separation 
between the shot and the sonobuoy, and V is the sound speed. 



In general, the separation X can be related to the depth H 



and the length of the reflected ray path by use of the Pythagorean 
Theorem; i .e . , 



where R is the distance from the shot point to the image point 0 1 . 

As successive shots are fired while steadily increasing 
the separation, the record of the direct travel times results in a straight 
line on the T/X plot shown in the upper half of Figure 2 . From Equation 
(5) it can be seen that the slope of this line is the reciprocal of the 
water sound speed. The record of the arrival of the reflected rays form 
a hyperbola asymptotic to the direct travel time line at large values 




( 6 ) 
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of X. Equation (6) can be combined with the relation 




where T is the travel time of the reflected ray and V is the sound speed 
to obtain the equation of a hyperbola 



shown that an additional term is introduced into the above equation. 
The equation for the reflected travel time for such a model is 



where 9 is the angle of slope of the bottom. The sign of the additional 
term is negative if the bottom dips in the direction opposite to that 
shown in Figure 3 . 



sub-bottom reflector is considered, it will be realized that the record 
of reflection times from that reflector will also form a hyperbolic trace 
on the T/X graph in accordance with Equation (7) , with V equal to the 
mean sound speed between the surface and the reflector and H equal 
to the depth to the reflector. The mean sound speed V can be calculated, 
but it provides little information for determining the sound speed-depth 




(7) 



For a bottom with a slope shown in Figure 3 , it can be 



12 2 1/2 
T = — (X + 4H + 4HX sin 9) 7 



( 8 ) 



If a mean sound speed between the sea surface and any 
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Figure 3 . Ray paths and T/X curves for rays 

reflected from a plane dipping surface 
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relationship. The sound speed that is required for determining this 

relationship, however, is the interval sound speed of the layer defined 

by an upper and lower reflector. 

2 . Determination of Interval Sound Speed 

The determination of the interval sound speeds by reflection 

measurements is an old geophysical problem (Clay and Rona , 1965) and 

techniques for reducing the data go back a number of years (Green, 1938) . 

The method of determining the interval sound speeds between successive 

reflecting layers from wide-angle reflection data is presently quite well 

known and has been described extensively in the literature (Dix, 1955). 

Its use at sea has been discussed by Houtz and Ewing (1963), Knox 

(1965), Clay and Rona (1965), and Le Pichon, Ewing, and Houtz (1968), 

among others. The method used in this paper follows the technique 

introduced by Le Pichon, Ewing, and Houtz (1968) and is the method 

used in the computer program of Appendix A. 

The water layer is assumed to have a known sound speed- 

depth profile, with the horizontal sound speed at the surface V. dif — 

h 

ferent from that of the sounding speed (the sound speed integrated 
over the depth H of the water layer) . 

It is further assumed that the dip of the bottom is known. 

The dip of any sub-bottom layer will be considered here to be the angle 
between the upper and lower reflectors defining the layer. Dip is con- 
sidered positive if the distance between the two reflectors increases 
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in the direction of separation between source and receiver. The dip 



of the water layer illustrated in Figure 3, for example, is positive. 

Assuming such a dip and using a mean water sound speed 
equal to the sounding speed V , Equation (8) can be rewritten as 



t 2 = 



2 2 

(X + 4H + 4HXsin9) , 



(9) 



V, 



and Equation (1) becomes 



T = 
o 




( 10 ) 



Since the direct travel times are a function of the hori- 
zontal sound speed V, , Equation (5) becomes 

h 



D = 




( 11 ) 



Substituting Equation (10) and (11) into Equation (9) and 
rearranging terms, an equation relating the reflected travel time to the 
direct travel time is 



T 2 = T 2 + (— t 1 — ) 2 D 2 + 2T (— 7 1 — ) D sin 9. (12) 

o V 1 o V 

If is replaced by the mean vertical sound speed of the material con- 
tained between any two reflectors, called the interval sound speed 
V , Equation (12) will relate the reflected travel time within the 
layer to the direct travel time at the surface corresponding to the 
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reflected time. The subscript i+1 indicates the sound speed for the 
i-th sub-bottom layer. The wide-angle reflection records, however, 
relate total reflected travel time to total direct travel time. 

The technique used to solve for the interval sound speed 
involved reducing the travel times for effects of the upper layers in 
order that the reduced values have the form of Equation (12) . An 
iterative method is then used to reduce the effect of dip by eliminating 
the curvature term 0 of Equation (12); i.e. , 



0 





) sin 9. 



(13) 



Finally, a least-squares line is fitted to the reduced 
2 2 2 

(T - T )/D data and the coefficient 
o 



is determined . 
determined . 



V. 



V. 



i+1 



From this, the interval sound speed is readily 



a . First Layer 

The water between the surface of the ocean and the 
sea floor is considered here to be the first layer. It is assumed that 
the dip of the sea floor 9^ is known and the sounding speed can be 
determined with good accuracy. 

Knowing and 9^ and knowing the reflected travel 
times as a function of the direct travel times from the T/X data, the 
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solution of the first layer involves computing the thickness of the water 

layer H and the horizontal sound speed at the surface V . 

1 h 

If data are obtained close to T , a fourth -order 

o 

least-squares polynomial is fitted to the T/X data, 

2 3 4 

T = a + a.D + a„D + a 0 D + a .D , (14) 

o 1 2 3 4 

where the a's are coefficients of the polynomial. At reflected travel 

time, T , the direct travel time D is zero; the reflected travel time, 
o 

also being the minimum travel time for the reflected wave, is equal 

to the coefficient a ; i.e. , 

o 

T = a . (15) 

o o 

The fourth-order curve is then extrapolated, obtaining the minimum 

reflection time T . 

o 

Care must be used with this method because the 

least-squares fitting technique does not constrain the fourth-order 

curve beyond the limits of the data . A more reliable but less accurate 

2 2 

method is to fit the T /X data to a linear least-squares line of the form 

T 2 =b +b 0 D 2 . (16) 

o 2 

2 

Again it can be shown that the coefficient b is equal to T . Thus, 

o o 

T is obtained by extrapolating a straight line instead of a curve as 
o 

was done in the case of the fourth-order polynomial. 
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The computer program used to reduce the T/X data 

decides which of the above two methods will be employed to solve for 

T . Three-tenths of a second direct travel time was considered the 
o 

maximum allowable distance from which the curve could be extrapolated 

using the fourth-order polynomial. 

Having computed T q and knowing V , the depth of the 

layer is computed using Equation (1). 

The curvature term in Equation (12) is small compared 

to the other two terms and, therefore, little error is introduced by 

assuming V, /V. equal to unity. This assumption is used to define a 
h 1 

correction term to be applied to Equation (12); i.e. , 

* = T 2 + 2T (1) sin 0 (17) 

o o 

The reduced travel times T are then computed for each data point using 

r 

the relation 



2 7 

T = T - * 
r 



(18) 



Combining Equations (12) and (18) , the expression for T^ becomes 



Tr 2 = e>->V. 



(19) 
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where T is a reduced travel time with the effects of both the dip and T 
r o 

eliminated. The horizontal water sound speed V is then found by fitting 

h 

2 2 

T and D to a least-squares first-order line of the form 

r 

T r 2 = C;i D 2 , (20) 

and solving for the coefficient c^. Knowing V^, Equation (19) can be 

solved for V. ; i . e . , 
h 

v h =v i^T • < 21 > 

The separation X is then computed using Equation (5), 



X = DV l (22) 

h 

b . Second Layer 

The second layer is, in fact, the first layer of sedi- 
ments for which an interval sound speed is desired. The record of 
reflected travel times from the lower interface of the second layer T 
forms a hyperbola similar to that of the first layer. In the first layer 




puted for each point. In the second layer, the interval sound speed 
is unknown, but the reduced direct travel time can be deduced by reducing 
the effects of the first layer on the travel times of the second layer. 
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The thickness of the second layer H ^ is computed 

using assumed values for the interval sound speed and dip 9^ of 

a a 

the layer. These values are referred to as the trial solution for the 
second layer. The trial solution is based on the apparent dips read on 
the normal reflection profiles, converted to true dips using trial sound 
speeds. The exact value of these trial sound speeds are not critical, 
so long as they are consistent with the dips entered. 

In Figure 4, an assumed ray path is drawn as the 
dashed line SBC. This ray is perpendicular to the second reflector; i.e., 
to the lower interface of the second layer. From Snell's Law 

sin 0 = (— sin 9 ) . (23) 

V 2 Z 

The distance the ray travels in the first layer SB is 



therefore 



SB 




cos 0 



(24) 



and the time T^ 



in the first layer 



T r 



SB 

V, 



where and are both known. 



(25) 
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Figure 4. T q reflection from the second layer 
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To obtain the reduced minimum travel time T for the 

°2 

r 

second layer alone, the minimum measured time to the second reflector 

T (obtained fitting the data to a least-squares polynomial) must be 
°2 

reduced by the time the sound traveled in the first layer T^, 



T = T - T. . 

°2 °2 1 
r 



(2 6 ) 



The length BC is obtained using the trial sound speed 



V 2 ’ 
a 



BC = V 0 T 

2 °o 
a 2 



(27) 



The length AB is the distance the direct ray will travel 
corresponding to the distance SB the reflected ray travels in the first 
layer. AB is computed from the expression 



AB = H x tan 0 , (28) 

and the depth of the second layer below the shot point S at reflected 

travel time T is 
o 

H 2 = BC + AB sin 9 2 . (29) 

H 2 is a first approximation for the thickness of the second layer. 
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can be written as 



Since was computed for the first layer, Equation (14) 



T = d Q + dj X + d 2 X 2 + d 3 X 3 + d 4 X 4 , (30) 

where 

X = DV U , 
h 

and 



a. 

1 




Differentiating with respect to X, the horizontal 
separation, Equation (30) becomes 



dT 

dX 



2d 2 X 



+ 3d 3 X 



4d 4 X' 



( 31 ) 



It can be shown (Clay and Rona, 1965) that the angle of emergence 
at the sea surface j3+9, is related to the derivative by 



dT 

dX 



sin {0 + 9j) 



V 



1 



(32) 



The data from the T/X graph is fitted to a fourth-order 
least-square polynomial and the coefficients of Equation (30) determined. 
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1 



The angle (0+9) is computed for each data point by combining Equations 
(31) and (32) to give 



sin (0 + 0 ) = V (d + 2d 2 X + 3d 3 X 2 + 4d 4 X 3 ) . (33) 

Knowing the dip of the first layer 9 , 0 is easily obtained. 

From 0, and using the sounding speed V , the travel 
time of the second reflected ray within the first layer is computed. This 
is then subtracted from the total reflection time to obtain the reduced 

i 

reflected travel time T 2 , the time the reflected sound from the second 

r 

reflector travels within the second layer. Similarly, a reduced direct 

travel time D is obtained and from Equation (5) , the reduced separation 
r 

X . Thus, the effects of the first layer are eliminated for each point and 
the problem reduced to a single layer case. 

The effect of dip is removed in a way similar to that 
used for the first layer. Equation (17) becomes 



where 





) sin © 2 , 
a 



(34) 



X = V D , . (35) 

r hr 

and the subscript 'a' indicates the assumed values of the trial solution 
of V 2 and 9 . The reduced times corresponding to the flat layer case T 
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are obtained from 



T 2 = T 2 - * . (36) 

r r 

The reduced values are then used to compute the co- 
efficients of a straight line equation of the form 



2 2 
T, = k + k X , 
2 o 1 r 
r 



(37) 



where k should be zero because in Equation (36) eliminates T 
o o 

and 






(38) 



and V can be readily computed. The dip 9 is then corrected by 

La Lu 

a 

V 2 

tan 0 = — tan 9 . (39) 

2 a 

a 

The solutions for V and 0 then replace the original 

U La 

trial solutions and the process is iterated a second time. A final 
solution is obtained for V 2 , 9 ^, and H 2 and the computations move on 
to the next layer. 

Thus, in the first layer, knowing the sounding speed 
and the dip 0^, and having the record of reflected travel times T^ 
verses direct travel times D^, the value of the water depth and 
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horizontal sound speed were computed. In the second layer, the 
problem was reduced to finding the interval sound speed V and the 
thickness by eliminating the effects of the first layer and then re- 
ducing the effect of dip using a trial solution. In each successive layer, 
the problem is first reduced to a single layer case by eliminating the 
effects of all upper layers using the computed values of thicknesses 
and interval sound speeds of those layers, and then the effect of dip is 
subtracted using an iterative process and a trial solution. 

This method is both fast and accurate. The time to 
reduce the data from a typical station is usually less than 20 seconds 
on the IBM 360/67 digital computer and the theoretical accuracy obtained 
is better than 1 part in 10,000 if the proper dips are entered (Le Pichon, 
Ewing, and Houtz, 1968). The primary advantage of this method is that 
there is no requirement that the layers be parallel, but rather, they may 
be dipping with respect to one another. 

C . THE COMPUTER PROGRAM 

The HONDO I program used to reduce the wide-angle reflections 
is listed in Appendix C. It was adapted from a similar program used at 
Lamont-Doherty Geological Observatory (LGO) of Columbia University 
(Houtz, personal communication) and follows the mathematical technique 
just presented of solving for thickness and interval sound speed of 
successive sub-bottom layers. Certain modifications were necessary 
to make the program compatible with the wide-angle reflection system 
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developed and for use in conjunction with the IBM 360/67 digital 
computer at the Naval Postgraduate School (NPS) using FORTRAN IV 
language. A thorough description of the routine followed by the program 
is presented in Appendix B. 

Sample data were included with the LGO program and modified 
accordingly (Appendix D) to be used with the HONDO I program. The 
corresponding results (Appendix E) were compared with those obtained 
using the IBM 1130 digital computer at LGO and the original form of the 
program, SLOW I. No significant difference was observed between the 
two solutions. The results from both programs are given in Table I. 

1 . Data Reduction 

Instructions for entering the data on IBM cards are included 
on the first page of the HONDO I program list. A more delicate part of 
the data reduction consists of selecting and tracing the major reflectors 
from the precision seismic graphic record made during the field profiling. 
Here special care must be taken to avoid multiples that can occur and to 
use only those traces representing the events from true lithological re- 
flectors . The difference between true reflectors and multiples is shown 
in Figure 5. It was assumed that the interval sound speed increased 
with depth . The apparent reflection times from multiples at a given 
separation are greater than those of the corresponding true reflector 
because the multiple travels longer in the lower speed layer than the 
true reflector. Note also, the T/X trace of the multiple is nearly parallel 
to the true reflection from the same reflector. Results that are not 
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Table I 



Comparison of Solutions of SLOW I 
and HONDO I Computer Programs 





SLOW I 


HONDO I 


Layer 


(LGO) 


(NPS) 


FIRST LAYER: 

Reflection Time 


5.216 


5.216 


(seconds) 
Sound Speed 


1.490 


1.490 


(k m/ sec) 
Propagated Error 


+.0023 


+.0023 


(km/sec) 

Depth 


3885.9 


3885.9 


(meters) 

SECOND LAYER: 

Reflection Time 


5.805 


5.805 


(seconds) 
Sound Speed 


1 .797 


1 . 7861 


(km/sec) 
Propagated Error 


+ . 0608 


+.0597 


(km/ sec) 
Thickness 


529.22 


526.28 


(meters) 
THIRD LAYER: 

Reflection Time 


6.250 


6.250 


(seconds) 
Sound Speed 


1.770 


1.7875 


(km/sec) 
Propagated Error 


+.0951 


+.0952 


(km/sec) 

Thickness 


379.83 


398.16 


(meters) 











SEPARATION X 




Figure 5. Multiple reflections 
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consistent with the local geology present another indication that multiples 
were present and traced. After the true reflectors are determined, the 
wide-angle reflection record is digitized in accordance with the pro- 
cedure outlined in Appendix A. 

2 . Trial Solutions 

As discussed previously, the trial solutions are obtained 
from normal profiling records. The trial values of interval sound speed 
and dip for each successive layer are entered on IBM cards in accordance 
with the instructions listed in Appendix C. If the interfaces are hori- 
zontal, there are no corrections for dip in the computations and the values 
for interval sound speed are unimportant. If the interfaces are not hori- 
zontal, inaccurate estimates of the sound speeds lead to only small 
errors in the first iteration and become negligible by the second iteration. 
However, dips entered which are inconsistent with the trial interval sound 
speeds may result in large errors. The error caused by entering such 
inconsistent values was investigated and the results given in Table II. 

The true solutions were those values computed using the data as given 
in Appendix D. The values computed in Case A were obtained assuming 
a trial sound speed 0.1 km/sec greater than the given data for both the 
second and third layers. For Case B, the trial sound speeds were un- 
changed, but the dips of the second and third layers were entered as 
+1.0°. It can be seen from Table II that the error introduced by incorrect 
sound speeds is small, but entering dips inconsistent with the trial sound 
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speeds causes significant error. Furthermore, the computer program 
will not detect such errors but will compute values consistent with the 
dips entered. 

Table II 

Comparison of Results of 
Various Trial Solutions 

Computer Sound Speed (km/sec) Propagated Error (+km/sec) 





Second layer 


Third layer 


Second layer 


Third layer 


True Solution: 


1.7861 


1.7875 


0.06 


0.10 


Case A: 


1.7861 


1.7875 


0.06 


0.10 


Case B: 


1.5756 


1.3724 


0.07 


0.08 



The total time to compute the results (Appendix E) from the sample 
data (Appendix D) was 16.36 seconds on the IBM 360/67 digital computer. 
The accuracy of the results is considered to be better than 1 part in 10,000 
if the correct dips are entered. The propagated error gives no indication 
that the correct dips were entered or that the computed values represent 
the true interval sound speed. Only by subjectively analyzing the results 
can gross errors be realized. 
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IV. INSTRUMENTATION 



A. THE WIDE-ANGLE REFLECTION SYSTEM 

The system for recording wide-angle reflection data is shown 
schematically in Figure 6. The profile recorder, timer switch, sound 
source, amplifiers and filters are those components normally carried 
on board oceanographic survey vessels for collecting normal-reflection 
seismic profiles. For collecting wide-angle seismic reflections, a 
sonobuoy system was developed to be used in conjunction with the on- 
board components. The sonobuoy system included a radio sonobuoy, an 
antenna, and a radio receiver. The system was developed with emphasis 
on simplicity and low cost; however, the two were not always synonymous. 

1 . The Sonobuoy 

The sonobuoy, launched while the ship was underway, per- 
formed as a self-contained station for detecting the sound pulses reaching 
that point in the ocean by various paths as shown in Figure 1. The 
pressure pulses were received by the hydrophone suspended below the 
sonobuoy by a thin wire. The pulses were converted to electrical signals 
and transmitted to the ship by the VHF radio transmitter contained in the 
watertight compartment of the sonobuoy. The sonobuoy transmitted on a 
pre-set frequency indicated by the channel number painted on the case. 
Table III lists the 31 sonobuoy channels available and the corresponding 
frequencies . The transmitter and acoustical amplifiers are powered by 
the sea-water-activated batteries located in the base of the housing 
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Figure 6. System for collecting wide-angle reflections 
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sound 

source 



and become active within seconds after the sonobuoy contacts the water. 
The sonobuoys used were of the AN/SSQ-23A, AN/SSQ-41, and 
AN/SSQ-5 7 types. Figure 7 is a drawing of an AN/SSQ-57 sonobuoy 
as it would appear in the water showing the antenna and hydrophone 
package . The characteristics of the various types of sonobuoys used 
are listed in Table IV. 



Table III 



Sonobuoy Frequencies 



annel 


Frequency 


Channel 


Frequency 




MHz 




MHz 


1 


162.25 


17 


162 .625 


2 


163.00 


18 


163.375 


3 


163.75 


19 


164.125 


4 


164.50 


20 


164 .875 


5 


165.25 


21 


165.625 


6 


166.00 


22 


166.375 


7 


166.75 


23 


167.125 


8 


167.50 


24 


167.875 


9 


168.25 


25 


168.625 


10 


169.00 


26* 


169.375 


11 


169.75 


2 7* 


170.125 


12* 


170.50 


2 8* 


170.875 


13 


171.25 


29 


171.625 


14 


172.00 


30* 


172 .375 


15 


172.75 


31* 


173.125 


16* 


173.50 






* crystals available for radio 


receiver 
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Figure 7. AN/SSQ-57 Sonobuoy 
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Table IV 





Sonobuoy Characteristics 

AN/SSQ-23A AN/SSQ-41 


AN/SSQ-57 


Weight, 

Kg 


8.11 


9.02 


9.02 


Channels 


1 - 16 


1 - 31 


1 - 31 


Time Selection, 
Hours 


1 


1 or 3 


1 or 8 


Depth of Hydrophone 
Meters 


, 18.2 


18.2 or 90.1 


28.9 


Transmitter Power, 
Watts 


1/2 


1 


1 



Each sonobuoy has a water-soluble scuttling plug located 
above the center of the housing. When the plug dissolves, water enters 
the watertight compartment and the buoy sinks. The floating period of 
the sonobuoy is a function of water temperature, with the plug dissolving 
more rapidly in warmer water. In no case, however, is the floating 
period less than 8 hr nor more than 20 hr in the case of the AN/SSQ-57 
type (Naval Air Systems Command, 1967). 

The plug was sealed with waterproof tape when a study was 
conducted investigating modifications to be made to the sonobuoy for 
recovering and re-using it on subsequent stations. The study revealed, 
however, that the rubber gaskets used to seal the water tight compartment 
were not intended for long use and that small amounts of sea water entered 
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the compartment after only 3 or 4 hr of immersion. This was learned after 



a considerable amount of sea water was found inside one particular 
sonobuoy after it was in the water for about 12 hr. On another occasion, 
a sonobuoy operated successfully for 3 hr in the water before it was 
recovered. The salt water batteries were removed immediately to prevent 
any electrolysis or heating damage and the buoy was washed in a mild 
solution of soap and water. Later, however, it was found to be inoperable 
due to corrosion of the electronic components inside the watertight case. 

In each case, the scuttling plugs were sealed. Considering this, and 
also considering the problem of trying to recover 18 to 90 m of hydrophone 
wire, which had a strong tendency to kink, it was decided that the cost 
of using one sonobuoy per station was far less than the complications 
which would incur and the time involved to recover and re-use a sonobuoy. 

The sonobuoys were launched by hand from the fantail of a 
ship. A screwdriver was used to bend the ends of the retaining ring on 
the bottom of the sonobuoy and release the bottom plate (Fig. 8). In 
some cases, the small screw used to secure the ring was removed and 
the entire retaining ring removed. Once the ring was removed, the 
bottom plate and hydrophone package were free to drop out of the sonobuoy, 
but were held in place until the sonobuoy was launched. The bottom 
plate serves as a weight to pull an aluminum protective housing away 
from the hydrophone package once the sonobuoy is in the water. For this 
reason, the clips (Fig. 9), which secure the plate to the housing, are not 
removed . 
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Figure 8. Releasing the bottom plate of the 
AN/SSQ-5 7 sonobuoy 
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Figure 9. Bottom plate retaining clips 
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The next step was to release a second retainer ring 
located around and near the top of the sonobuoy. Considerable caution 
was warranted here as this ring allowed the rotochute assembly to fall 
free and pull a small white lanyard releasing the spring-loaded transmitting 
antenna. The antenna erects with sufficient force to cause injury. Once 
it was safe to do so, the retaining ring was released by lifting a small 
U-shaped pin (Fig. 10) and the rotochute assembly gently removed. A 
slight tug given on the lanyard tripped the antenna retainer clip allowing 
the antenna to extend to its full operating position shown in Figure 7. 

The sonobuoy was then ready to be launched. 

2 . Receiver System 

A radio receiver and antenna system was used to receive 
and amplify the VHF signals transmitted by the sonobuoy. The radio 
used for this was a Japanese-made Realistic Patrolman pro-2 solid state 
VHF receiver. The decision to use this radio was based mainly on its 
availability at the time the study was started and the fact that its frequency 
band covered that of the entire sonobuoy range. Also of importance, and 
in keeping with the idea of simplicity and economy, the pro-2 was small, 
easy to carry, simple to operate, and relatively inexpensive ($99.00). 

The pro-2 radio receiver is a self-contained unit operated 
on normal 110-V AC. An internal speaker was used to monitor the signals. 
The signal was tuned by manual control for maximum volume. Crystal 
operation was available, but did not appear to increase the quality of 
the signal. The crystal did, however, provide a means for tuning the 
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Figure 10. Removing the rotochute assembly 
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receiver without a signal present. With the proper crystal in place, the 



radio was simply tuned for maximum background noise. Crystals for the 
seven sonobuoy channels indicated in Table III were purchased from 
Radio Shack Warehouse, 2615 West 7th Street, Fort Worth, Texas. 

The pro-2 receiver had two locations where the signal was 
available for recording. The signal taken from the speaker jack on the 
front panel of the cabinet was found generally to be highly distorted by 
the audio amplifier section of the receiver. The tape-out jack on the 
back of the cabinet provided a signal of rather low voltage, but when 
amplified it was found to be relatively undistorted . The results were 
further verified by using an oscilloscope to observe the wave form of 
the signal from each jack. A pure tone was transmitted from a signal 
generator and received by the radio. The signal at the tape-out jack 
was clear and undistorted while that of the speaker jack was either 
distorted at high volume settings or was lost in the amplifier noise at 
low settings . 

An 11-element Winegard Yagi antenna (Fig. 11) was used 
with the pro-2 radio receiver for receiving the signal transmitted by the 
sonobuoy. The Yagi was purchased when the ground plane dipole antenna 
(salvaged from a used sonobuoy) failed to provide sufficient power gain 
for receiving the signal from the sonobuoy except at very small separation 
distances . 
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Figure 11. Yagi 1 1 -element antenna 






49 



The antenna patterns for each antenna were recorded over a 



range of frequencies covering the sonobuoy range (Fig. 12). It can be 
seen from this graph that the Yagi was quite frequency-dependent and 
at lower frequencies did not produce better results than were obtained 
using the dipole antenna. Figure 13 shows the antenna pattern for the 
Yagi at various frequencies. The dipole antenna pattern was omni- 
directional over the entire sonobuoy frequency range. It was decided, 
therefore, that only higher frequency sonobuoys would be used and 
special attention would be given to insure that the Yagi antenna was 
pointed in the direction of the sonobuoy. 

A 300/50-ohm impedance matching transformer was attached 
at the antenna to couple the 300-ohm antenna to the 50-ohm shielded 
cable. All external connections were covered with alternate layers of 
tape and electrical coating to prevent salt spray from damaging the system. 
A stainless steel pipe, 7-m long and 3.75 cm in diameter, was used for 
an antenna support. 

The entire wide-angle reflection system weighed less than 
10 kg per sonobuoy and was quite transportable; the whole system can 
be carried in the backseat of an automobile. It took only a few minutes 
to set up the equipment aboard ship. The system was simple to use, 
requiring little more than to launch a sonobuoy and to tune the receiver. 
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POWER GAIN RELATIVE TO 10 db 




Figure 12. Power gain vs. frequency curve 
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Figure 13. Antenna patterns for Yagi antenna at 

various frequencies (relative to 175 MHz) 
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V. CONCLUSIONS 



It was found from this study that the Realistic pro-2 radio receiver, 
when used with a Winegard Yagi antenna, could adequately receive the 
VHF signal transmitted by a radio sonobuoy and at sufficient distance to 
conduct wide-angle seismic reflection profiling. It was further found 
that the signal received with this system was distorted by the audio 
amplifier section of the radio. An undistorted signal was available at the 
tape-out connection on the back panel of the pro-2 receiver cabinet, by- 
passing the audio amplifier section. 

It was concluded after becoming familiar with three types of radio 
sonobuoys that, when used normally, they were quite simple to operate; 
but any modifications that were possible required a considerable amount 
of preparation in order to use the sonobuoy for more than one station. 

Magnetic tape recordings have been made of the signal received 
using the sonobuoy system and the sample set of data from LGO demon- 
strated that the HONDO I computer program adapted for use with the 
system was able to reduce wide-angle reflections and obtain the interval 
sound speeds and corresponding thicknesses of sub-bottom layers. 

It is concluded, therefore, that the system developed can be used 
to collect wide-angle seismic reflection data and that the HONDO I 
program can be used to reduce the data. 
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APPENDIX A 



DIGITIZING THE WIDE-ANGLE REFLECTION RECORD 

If T is marked on the record (the instant when the sonobuoy and 
o 

sound source are together), this point becomes the origin. Otherwise, 
the direct travel time line, called the direct line, is extended to intercept 
the X-axis, thus defining the origin. A plastic overlay is placed over 
the record so as to include the origin and all direct and reflected traces 
within the millimeter-size grid etched on it. 

Figure 14 shows a sample record with a direct line and one reflected 
trace. The origin of the grid does not coincide with the origin of the 
record, but rather the entire data is included under the grid. Below the 
record are sample values read from the record, including the two points 
A and B used to set the time scale. Points A and B are selected anywhere 
along the T-coordinate and are a measure of the distance in millimeters 
from the origin of the grid in the T-direction only. The time DDD rep- 
resents the known time interval between the two points. 

The values from the direct line are recorded in eight-digit integer 
form. For example, a point located at 12 7.9 mm from the grid origin on 
the X -scale and 148.9 mm on the T- scale is recorded as 12791489. The 
X-scale is the vertical scale and represents the first four digits , while 
the T-scale is horizontal and represents the last four. From two to six 
points can be used to represent the direct line. The reflection traces 
are digitized similarly, except that as many as 90 points can be used. 

The program can handle up to 10 reflecting horizons. 
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TIME-AXIS (mm) 



Direct line: 



11101049 12791489 12931056 13781248 



Reflected line: 



11041373 11291382 11531395 11811414 12801438 
12281458 12501480 12611493 12811060 13001095 
13111120 



Point A: 



110.0 



Point B: 



140.0 



DDD: 



3 seconds 



Figure 14. Digitizing the data 
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At times, the sound pulse is detected at the sonobuoy after the 
recorder completes a cycle. This occurs when the separation between 
the sonobuoy and the sound source becomes large enough so that the 
travel time for a particular sound wave is greater than the sweep rate 
of the recorder. An example of this is shown in Figure 14 where the 
trace goes off the bottom of the record and continues on the top during 
the next cycle. In this case, the shift is taken care of v/ithin the 
program . 
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APPENDIX B 



ROUTINE FOLLOWED BY HONDO I COMPUTER PROGRAM 

The computations begin in the MAIN routine by reading from the 
first card in the data deck the station identification, ADEN1 and ADEN 2 , 
and the number of layers or reflections observed, NREF, for the particular 
station. Subroutine TRIALS is then called to enter the valves for the trial 
solutions which include trial sound speeds V from the second card in the 
data deck and the trial dips W from the third card. Only NREF values are 
taken from each card and blank spaces are read as zero. As an example, 
for the sample station DIANNE (Appendix D) , the card representing the 
trial dips is blank; therefore, the trial dip for each of the three layers 
is zero. The technique of reading zeros from blank cards is used through- 
out the program to control the exits from particular subroutines. For this 
reason, it is extremely important that a blank card follow the cards con- 
taining the digitized data pairs for the direct trace and for each reflected 
trace, and that an additional blank card follow the data from the last 
layer of the last station investigated. This additional blank card causes 
the program to exit the computer and completes the job. The blank card 
after each set of data points tells the computer that there are no more 
points for that particular trace. 

The fourth card in the data deck is also read by TRIALS. The card 
contains the value of the computed sounding speed W for the water layer. 
This value will replace the value read in as a trial speed for this layer. 
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Thus, the first number in the array labeled V is the actual value used for 



the sounding speed VV of the first, or water, layer. Also, the dip read 
for the first layer is not, in fact, a trial dip; but rather it represents the 
actual value determined for the sea floor. For this reason, the values 
entered for the first layer of sounding speed W and dip W must be those 
values computed from observations made at the time the data were col- 
lected . 

Values for V and W are transferred back to MAIN through COMMON. 

The program is returned to MAIN and the main loop for NREF layers 
begins by setting the integer L equal to 1 for the first layer. MAIN then 
calls both subroutines READS and INTERP in succession. 

For the first layer, READS enters the values from the fifth data card 
representing the constants which appear in Figure 14. They include two 
arbitrary points on the time scale A and B, the time in seconds represented 
by the distance between them (DDD) , the seconds required for each sweep 
of the recorder (SPS) , and finally, DELAY, which is a provision to include 
any time delay which might occur between the shot and the start of the 
recorder sweep. READS uses A, B, and DDD to scale the time axis in 
seconds. SPS is used to add the appropriate number of seconds to each 
point of the trace if the trace is continued at the top of the record during 
a scale shift. 

Once the time axis is scaled, READS then enters the digitized data 
points from the direct travel time trace. These points are converted to DX 
in millimeters and DY in seconds from the origin. 
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Next, READS enters the digitized points from the first reflected 



time trace and converts these points to RX and RY (in millimeters and 
seconds, respectively, from the origin). 

Subroutine COUNT is called from subroutine READS to count the 
number of non-zero data points found on each card. The total number of 
points for any reflected time trace is summed on II and carried via 
COMMON to other locations in the program. ND is the number of points 
used to compute the direct time line. 

All values entered in READS are transferred to MAIN and then to 
INTERP through COMMON. 

INTERP interprets RX in millimeters for DD in seconds of direct 
travel time corresponding to each reflected travel time RY. This is 
accomplished by computing the slope of the direct time line, DELXY, 
and multiplying it by the millimeters of separation between RX and the 
origin DX1 as shown in Figure 15. RX is then reset equal to DD at each 
point on the reflected trace. Thus, RX becomes the travel time corresponding 
to each reflected travel time RY. RX and RY are transferred to MAIN 
through COMMON as U and S respectively. In MAIN, the direct travel 
time corresponding to the first reflection time recorded is compared with 
the minimum allowable separation for considering a fourth-order polynomial 
fit to the data, DMIN . DMIN is equal to 0 . 3 sec. If the direct time is 
less than DMIN, subroutine LSFIT is called and the data are used to 
compute a fourth-order polynomial using a least-squares technique. If 
the data are not obtained at or near vertical incidence and the minimum 
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X 




DD = DY1 + (RX - DX1 ) * DELXY 



Figure 15. Relating direct travel time to 

corresponding reflected travel time 
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reflection time is greater than DMIN , the data are squared and a linear 



least-square polynomial is computed in subroutine LSFIT. In either case, 
the data are extrapolated to T q using the resulting polynomial equation,, 
the minimum reflection time is obtained, and the depth is computed using 
the sounding speed. 

Subroutine LSFIT also computes the theoretical reflection time 

corresponding to the measured direct time by solving for the equation of 

2 2 

a line fitted to the T /X data-. 

Back in the MAIN routine, the corrective term CON to be applied 
to the travel time is computed and subtracted from the total time for 
each point. Thus, the reduced times R and D are squared and LSFIT is 
called to fit the values to a linear least-squares line and compute the 
horizontal sound speed from the coefficient X (2) . 

This completes the computations for the first layer, having obtained 
the depth and horizontal sound speed. These results are printed and the 
program returns setting L equal to 2 . 

In successive layers, the program begins by calling the subroutine 
READS. The digitized data for the particular layer is read and converted 
to millimeters of separation and seconds of travel time from the origin. 
Subroutine INTERP is called, interpreting millimeters of separation for 
seconds of direct travel time as was done in the case of the first layer. 
Now, however, the horizontal sound speed is known and the separation 
can be converted to distance in kilometers by an equation equivalent to 
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Equation (5). The values of reflected times and corresponding values of 
horizontal separation are returned to MAIN through COMMON as S and U 
respectively. 

In MAIN, the horizontal separation for the first point is compared 
with the minimum allowable separation for considering a fourth-order 
polynomial, DMIN. In this case, DMIN is 0.45 km, obtained from the 
product of an assumed horizontal sound speed of 1 . 5 km/sec and a 
minimum allowable separation. time of 0.3 sec. If the first point is less 
than DMIN, the data are used in the subroutine LSFIT to compute the co- 
efficients of a fourth-order least-squares polynomial. If the first point 
is greater than DMIN, the values are squared and used to compute a 
linear lease-squares equation, again in subroutine LSFIT. As was done 
for the first layer, the data are extrapolated to obtain the minimum re- 
flection time for the particular layer being considered. 

The derivative of the fourth-order polynomial with respect to the 
horizontal separation is used to compute the angle of emergence 
ZM {0 + 9 in Fig. 4) using an equation equivalent to Equation (33) . If 
the minimum reflection time was previously computed using a linear fit 
to the data, subroutine LSFIT is called a second time to compute the 
coefficients of the fourth-order polynomial fit. ZM is a function of the 
slope (dT/dX) of the fourth-order curve. The curve is not restrained 
beyond the first and last points and its slope may become erratic there. 
For this reason, the angle of emergence is not computed for the first 
and last points . 



62 



Subroutine TUPPER is called and the thickness of the layer being 



considered is computed using the trial solution of interval sound speed 
and dip. The technique follows that given in the previous section for 
the first layer, using an equation equivalent to Equation (29) . 

Subroutine TUPPER uses the emergence angles to compute travel 
times TT in the layers above the layer for which a solution is sought. 

TT is computed using the previously-solved values of interval sound 
speed and dip for the upper layers. The corresponding values of hori- 
zontal separation are found in a similar way. The travel times and 
separations are transferred to MAIN as TR and DR respectively. 

MAIN reduces the measured values of reflection time S and 
horizontal separation U to obtain D, the reduced value of horizontal 
separation, and R, the reduced reflected travel time. D is also adjusted 
for the sum of the dips of the upper layers by dividing by C, the cosine 
of the sum. 

The reduced D and R corresponds to a single layer case as shown 
in Figure 3 . 

The travel time for each point is finally reduced to that of a flat- 
layer case by subtracting equations equivalent to Equations (10) and (15). 
The value R in the computer program is now the reduced time squared for 
the flat-layer case. The value D is also squared and the subroutine 
LSFIT is called to fit a line to the reduced values. The second coef- 
ficient X(2) of the linear equation obtained from the subroutine LSFIT is 
then used to compute the interval sound speed. 
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The dip of the layer W is corrected using the relation 



V 

tan W = rr tan W , 

V Q 
a 

where the subscript 'a' indicates the values for the assumed trial 
solution . 

This solution then replaces the original trial solution and the 
computations go back to computing a new thickness and start the second 
iteration. Upon completion of the second iteration, the solutions for 
the layer are printed and the routine starts a new layer. 

The results are tabulated and printed for each layer in succession. 

A sample list of results for the theoretical station DIANNE are included 
in Appendix E . 

The first line gives the station name, layer number, and number 

of points used by the programs to compute the results. The number of 

points considered may be less than are digitized due to rejection of 

points occurring before T or those resulting in negative reduced times. 

o 

The first three columns of figures below the heading represent squared 
values of reduced separation and reduced reflected travel times cor- 
responding to a flat-layer case. For the first layer, the square of the 
separation listed under the column labeled X2 , are in seconds squared, 
whereas, for each of the following layers, they are in kilometers squared. 

The columns labeled T2 MEASURED and T2 COMPUTED are the 
reflection times squared corresponding to the separation in the first 
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column. The measured values are those values read from the reflection 

curves and reduced to the flat-layer case. The computed values are 

from a least-squares curve fitted to the reduced data. The fourth column 

is the difference between the measured and computed values of the 

reduced reflection times . 

2 

The variance cr of the data is computed from 



2 



cr 



l 

N-2 




<T 2 - T V 

c 



where N is the number of points, T the observed times, and T the 

c 

computed times. The standard deviation is the positive square root of 
the variance. 

If many reflection profiles were to be taken at the same location, 
the variance of the interval sound speed about the mean can be calculated. 
However, with only a single set of data, the variance cannot be calculated 
in this way. Instead, the techniques used to compute the variance, or 
propagated error, for each interval sound speed is based on the variance 
of the coefficients for each point compared to those coefficients obtained 
with the least-squares fit. 
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APPENDIX C. HONDO I COMPUTER PROGRAM 
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cc 
cc 
cc 
cc 
cc 
cc 
cc 
cc 
cc 
cc 
cc 

WAS PREPARED BY CC 

OF LAMONT-DOHERTY GEOLOGICAL OBSER-CC 
VATORY OF COLOMBIA UNIVERSITY. THE PRESENT PROGRAM CC 
WAS ADAPTED FOR USE ON THE IBM 360 BY S.K. EDLESON CC 
LCDR,USM AT THE NAVAL POSTGRADUATE SCHOOL, MONTEREY , C-^ 



HONDO I 

PROGRAM TO SOLVE FOR THE DEPTH, DIP AND INTERVAL 
SOUND SPEEDS WITHIN THE SUBSURFACE LAYERS. 

THIS IS ACCOMPLISHED BY SOLVING THE WIDE ANGLE RE- 
FLECTION PROFILE 
FLAT LAYER CASE. 

ORIGINAL PROGRAM 
XAVIER LE PICHON 



BY REDUCING EACH LAYER TO A SINGLE 



JULY, 1970. 



CARD 



THE 

ARE 



FOLLOWING 

REQUIRED: 



INPUTS 



COLUMNS 



1-8 STATION IDENTIFICATION 

14-15 NUMBER OF LAYERS (RIGHT ADJUSTED) 

1-10, ETC TRIAL SOUND SPEEDS FOR EACH LAYER IN 
KM/SEC. MORE THAN ONE CARD MAY BE 
USED , BUT IN NO CASE CAN THERE BE 
MORE THAN TEN LAYERS. 

1-10, ETC TRIAL DIPS (DEGREES PLUS OR MINUS) 



1-10 

1-6 

7-12 

13-18 

19-24 

25-30 

3-10, 
13-20, 
23-30 , 
33-40, 
43-50 , 
ETC. 



1-8 

11-15 



3-10, ETC 



COMPUTED SOUNDING SPEED IN KM/SEC 

ANY REFERENCE MARKING ON R-TIME SCALE. 
(ENTERED IN MM AND TENTHS, I. E. 102. 4MM 
ANOTHER MARK ON R-TIME SCALE. 

NUMBER OF SECONDS BETWEEN ABOVE TWO 
MARKI NGS 

NUMBER OF SECONDS PER CYCLE 

TIME DELAY BETWEEN RECORDER AND SHOT 

DIGITIZED DATA PAIRS READ OFF D-LINE. 
FIRST FOUR NUMERALS ARE MM 8 1/10 *10 
FROM X-TIMF SCALE AND SECOND FOUR THE 
SAME FROM THE R-TIME SCALE: I.E. IF A 

POINT IS (102.4, 237.1) IT IS ENTERED 
10242371. USE AS MANY CARDS AS NEEDED 
EACH POINT MUST BE GREATER THAN 100.0 



BLANK (DON’T FORGET BLANK 
ARE IMPORTANT) 



CARDS, THEY 



STATION IDENTIFICATION 
NUMBER OF LAYER THE FOLLOWING DATA 
REPRESENTS (RIGHT ADJUSTED). THIS CARD 
IS REQUIRED FOR EACH LAYER 

DIGITIZED DATA FROM ONE LAYERS R-LINE, 
ENTER SAME AS D-LINE DATA 

BLANK CARD REQUIRED BETWEEN EACH LAYER 
AND TWO BLANK CARDS ARE REQUIRED AFTER 
THE LAST LAYER OF THE LAST STATION. 



CC 

CC 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

ACC 

.cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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HONDO I 



MAIN ROUTINE 



IMPLICIT REAL*8( A-H,0-Z) 

DIMENSIONTOI 10) ,HH( 10) , V( 10) , W( 10) , SW ( 1 0) , XXXX( 2 5 ) , * 
ICW( 10) ,TW( 10) ,D( 90), R( 90),ZM( 90 ) , P ( 5 , 6 ) , X ( 5 ) 
DIMENSION U( 90), S< 90 ) , DX ( 2 5 ) , D Y ( 2 5 ) , T I ( 90 ) , DDL ( 1 0 ) 
COMMON T0,HH, V, W,TR, DR, SW,CW,TW,D, R, ZM, P,C , SVN, VH 
COMMON X , T I T ,U, S,XXXX, DY , DX0 , DY0 , YSCAL , SGN 

COMMON M , I NO , KO , I X X X X 
EQUI VALENCE(WW,DMI N) 

C READ IDENTIFICATION AND NUMBER LAYERS STATION 

1 READ! 5 , 1 10) ADEN1 , ADEN2 ,NREF 
110 FORMAT ( 2A4,2X, I 5) 

IF (NREF.EQ.O) GO TO 9999 

EXIT OUT OF ROUTINE IF NREF IS ZERO. REQUIRES LAST 
CARD BE BLANK. *****N0T E**** LAST TWO CARDS IN 

DATA DECK MUST BE BLANK. 

10 CALL TRIALS ( KW , ADEN1 , A DEN2 , NR E F , V V , PE RC E ) 

MAIN DO LOOP FOR NREF LAYERS 
READ IDENTIFICATION LAYER 
DO 12 L=1 , NREF 

401 CALL READS(L,ADEN1,ADEN2,IRFF,SPS) 

CALL INTERP ( PE RC E , A DEN L , ADEN2 , I REF , L , SPS ) 

FIT CURVE TO DATA TO OBTAIN T-ZERO AND LATER 
DIFFERENTIATE FOR ANGLE OF EMERGENCE 
IND IS AN INDEX USED IN LSFIT SUBROUTINE 
U IS NUMBER OF SECONDS OF DIRECT TRAVEL TIME CORRES- 
PONDING TO S SECONDS OF REFLECTED TRAVE TIME. 

402 IF(L-l) 930,931,930 

CHECK FOR DATA TO BE FITTED TO 4TH ORDER POLY: I.E. 

DATA MUST BE COLLECTED BEFORE .45 KM OR .3 SEC SEPARAT 
ION. UIl) IS SEPARATION FOR FIRST DATA POINT. 

930 DM I N = 0.45 
GO TO 932 

931 DM IN- 0.3 

932 IF (U(l).LE.DMIN) GO TO 3 
923 I ND=- 1 

DO 90 K= 1 , KO 

R AND D TAKE ON DIFFERENT VALUES CORRESPONDING TO 
REFLECTED AND DIRECT TRAVEL TIMES RESPECTFULLY 
R(K)=S(K)**2 
90 D(K)=U(K)**2 

CALL LSFIT (1,L) 

TOIL) -DSQRT (XII) ) 

I ND= 1 

I F ( L — 1 ) 2, 2,3 

IF FIRST LAYER COMPUTE THICKNESS AND VH BY T2/X2 

2 CON=2.*TO< 1)*DSIN (W(l)) 

HH ( 1 )=V( 1 )*T0< 1 ) /2. 

DO 200 1=1, KO 

BRING BACK TO FLAT LAYER CASE 
R ( I ) = S ( I )**2+C0N*U( I ) 

200 D( I )=U( I )**2 
I ND= 1 

CALL LSFIT (1,1) 

VH-DSQRT ( X ( 2 ) )*V(1 ) 

GO TO 32 

3 DO 92 K= 1 , KO 
R(K)=S(K) 

92 D(K)=U(K) 

CALL LSFIT (4,L) 

COMPUTE COEFFICIENTS OF THE DERIVITIVE FOR OBTAINING 
THE EMERGENCE ANGLE ZM 
X ( 2 ) =X ( 2 > - V ( 1) 

X(3)=X(3)*2.*V( 1) 

X( 4>=X(4)*3.*V( 1) 
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X( 5) =X( 5) *4.*V( 1 ) 

CHECK = D ( D-DMIN 
IF (CHECK. GT. 0.0) GO TO 920 

921 TO ( L ) = X ( 1 ) 

920 IF (L.LT.2) GO TO 2 

FIND ANGLE OF EMERGENCE ZM FOR EACH POINT EXCEPT 
FIRST AND LAST 

922 K0=K0— 1 
I ND= 1 

DO 30 1=2, KO 

ZM IS FIRST THE SINE OF THE ANGLE 

ZM( I ) =X(2)+X(3 >*U( I )+X(4)*U( I ) **2+X( 5) *U( I ) **3 

30 ZM ( I )=DARSIN(ZM( I ) ) 

ZM IS NOW THE ANGLE IN RADIANS. 

D ( 1 )=-l.D20 

T I ( 1 )=0. 

31 I ND=— I ND 

CON=2.*DSIN(W(L) ) / V ( L ) 

I K = 0 

DO 300 1=2, KO 

COMPUTE TRAVEL TIMES IN UPPER LAYERS 
CALL TUPPER (L ,ZM( I ) ,IK) 

I K = 1 

REDUCE SEPARATION U TO SEPARATION FOR L LAYER ONLY 
NOTE THAT D IS IN KM 
D ( I ) =U ( I ) -DR 

REDUCE REFLECTION TIMES S FOR TIMES ONLY WITHIN L 
R ( I ) =S ( I )-TR 
T I ( I ) = T I T 

NOTE TIT IS T-ZERO. REFLECTION TIME FOR LAYER L ONLY. 
CORRECT SEPARATIONS FOR DIPS OF UPPER LAYERS: 

RECALL C IS THE COSINE OF THE SUM OF THE DIPS. 

D( I )=D( I )/C 

ELIMINATE POINTS WITH NEGATIVE REDUCED D 
I F ( D ( I ) 1300,3601,3601 
REDUCE TO FLAT LAYER CASE 
3 601 R ( I ) = R ( I ) **2+C0N*D( I)*TI(J)-TI(I)**2 
D( I J=D ( I )**2 
300 CONTINUE 

SOLVE FOR SOUND SPEED 

KW= 1 FOR HORIZONTAL LAYERS: 

KW=2 FOR DIPPING LAYERS. 

GO TO ( 5,6) , KW 

5 I ND= 1 

6 CALL LSFIT ( 1,L) 

VV = V ( L ) 

C CHECK FOR IMAGINARY SPEED 

I F ( X ( 2 ) >3100,3101, 3101 
3101 V( L) =DSQRT( 1 . / X ( 2) ) 

C RESET ANGLE ACCORDING TO NEW SPEED 

TW(L )=TW( L ) * V ( L) /VV 
W( L) =D ATAN ( TW( L> ) 

C COMPUTE NEW THICKNESS HH. 

CALL TUPPER (L,ZM(2),0) 

C IF IND NEGATIVE RESTART PROCESS 

I F ( IND >31, 31, 32 

32 WW=W( L )*0. 572957795131D02 

WR ITE(6»104)V(L) , TO ( L ) ,SVN, WW 
IF(L.GT.l) GO TO 69 
DOL ( L)=V(L)*TO(L) *500.0 
WRITE(6, 105) DOL(L) 

105 FORMAT ( * S' DEPTH OF SEA FLOOR • , T 26 , ' = • , F8 . 2 , 

1T45, 'METERS' ) 

GO TO 12 

69 DOL ( L ) = V ( L ) * ( TO ( L ) -T0( L— 1 ) ) *500.0 
WR I T E ( 6,106) L , DOL ( L ) 

106 FORMAT ( ' ',' THICKNESS OF LAYER' ,1 3,T26, ' = ' ,F8.2, 

1 T45 , 'METERS' ) 

104 FORMAT ( ' ',' INTERVAL SP EE D ' ,T26 , • = ' , F8 . 4 , T45 , • KM/ S EC 
1'/' ',' T-ZERO REFLECTION TIME' ,T26, '=' ,F8.4,T45, 'SEC 

1'/' ',' PROPAGATED ERROR ', T26 ,' = +•/• + ' ,T28, '_' ,D15. 7, 

1T45, 'KM/SEC' / ' ',' DIP OF LAYER • ,T26 ,' = SF7.4, 
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1T45, 'DEGREES' ) 

12 CONTINUE 

WRITE(6,ill)VH 

111 FORMAT! '1',' SPEED OF SOUND AT WATER SURFACE 
1T37 , F7 . 3 , ' KM/SEC' ) 

GO TO 1 

3100 WRITE! 6,3200)1 

3200 FORMAT! • SPEED NEGATI VE , REMOVE LAYER', 13) 

9999 WRITE! 6,9000) 

9000 FORMAT! ' 1* ) 

STOP 

END 
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SUBROUTINE TRIALS 



C READ TRIAL SPEEDS AND CORRESPONDING DIPS IN DEGREES: 



SUBROUTINE TRIALS ( KW , ADEN I , ADEN2 , NR EF , VV , P ERC E ) 
IMPLICIT REAL*8( A-H,0-Z) 

DIME NS I ONTO! 10),HH(10),V!10),W(10) 

COMMON TO,HH,V,W 
10 READ! 5 , 107 ) ( V( I ) , 1=1 , NREF) 

READ ( 5, 107) (W( II , 1=1, NREF ) 

107 FORMAT ( 8F10. 5) 

KW=2 

WT=0. 

DO 6000 1=1, NREF 

6000 WT=WT+DABS( W( I ) ) 

SET K W= 1 IF LAYERS ARE HORIZONTAL: 

USED IN HONDO I . 

I F ( WT-l.D- 10)6001, 600 1,6002 

6001 KW= 1 

6002 WRITE! 6,103) ADEN1, ADEN2, NREF 

103 FORM AT ( 1H1 , ' STATION ' ,2A4, ' NUMBER LAYERS ',13/// 

WRI TE ( 6 , 600 ) 

600 FORMAT!//* TRIAL SPEEDS AND DIPS •/) 

WRITE! 6,500) 

500 FORMAT! 'O' ,RX, 'TRIAL* , 3X , ' TR I AL ' / • ', 

1T2, 'LAYER* ,3X, ' SPD. ' ,4X, 'DIP */' ', , 

IT 10, 'KM/SEC , 2 X , • DEGREES •/' + ' , 

1T2,5( '_') ,3X,6( ) ,2X,7('_')/) 

DO 555 1=1, NREF 
555 WRITE! 6,888) I , V < I ) , W ( I ) 

888 FORMAT!' ' ,I3,F11.4,F6.2) 

C CONVERT ANGLES TO RADIANS 

DO 1300 1=1, NREF 

1300 W< I )=W( I ) /O. 5729577951 31 D02 
C READ MAIN CONSTANTS 

READ! 5, 105) VV, PERCE 

105 FORMAT ( 2F 10. 5 ) 

C VV IS SOUNDING SPEED 

V! 1 ) =V V 

IF ( PERCE. GT.O. ID-10) GO TO 1105 
1100 PERCE =1.0 
1105 WRITE! 6, 106) VV 

106 FORMAT ( '0 1 , ’ SONIC SOUNDING SPEED =',F6.3,' KM/SEC') 
RETURN 

END 
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SUBROUTINE READS 



THIS SUBROUTINE READS DIGITIZED TABLE OUTPUT OF 
D-LINE AND R-CURVES. CONVERTS RY AND DY TO SECONDS 
RX AND RY ARF TRANSMITTED TO HONDO I BY COMMON AS 
U AND S RESPECTFULLY. 



SUBROUTINE RE ADS ( L , ADE N1 , ADE N2 , I RE F , SP S ) 

IMPLICIT REAL*8( A-H,0-Z) 

DIMENSION RX ( 90>,RY( 90 ) , DX ( 2 5 ) , DY ( 2 5 ) , I T ( 1 6 ) 

DIME NS I ONTO! 10),HH(10),V(10),W(10),SW(10), 

1CW<10) »TW< 10) ,D( 90), R< 90),ZM( 90 ) , P ( 5 , 6 ) , X ( 5 ) 
COMMON TO , HH , V , W , TR , DR , SW , CW , T W , D, R , ZM , P , C , SVN , VH 
COMMON X,TIT,RX,RY , DX , DY , DXO ,DYO , Y SC AL , SGN 
COMMON M, IND, I I ,ND 

IF THE FIRST LAYER, READ IN THE D-POINTS IN ORDER TO 
SCALE THE REFLECTED TIME SCALE 

IF(L-l) 1,1,3 
I 11=0 

RE AD ( 5 , 1002 ) A, B , ODD, SPS , DELAY 

1002 FORMAT (5F6.0) 

READ (5,889) (IT<I), 1=1, 12) 

889 FORMAT (6(2X,2I4)) 

I NDE X= 1 
IT ( 13) =0 
DXO= ITU) 

COMPUTE YSCAL EQUAL TO NUMBER OF SECONDS OF RECORD 
PER MILLIMETER OF MEASUREMENT OF THE REFLECTED 
TIME SCALE 
YSCAL=DDD/(A-B» 

DYO= I T { 2 ) 

20 CHECK = I T ( 3 ) - I T ( 1 ) 

IF (CHECK. GT. 0.0) GO TO 22 

21 SGN=- 1 . 

GO TO 5 

22 SGN= 1 « 

GO TO 5 

3 11 = 0 

I NDE X= 2 

READ ( 5,1003) AD EN 1 , AD EN2 , I REF , I R, DRE F 

1003 FORMAT ( 2 A4 » 2X ,2 15, F10.3) 

4 READ ( 5 , 1004) ( I T ( I ) , I = 1 , 16 ) 

1004 FORMAT ( 8 ( 2 X , 2 I 4 ) ) 

I F ( I T ( 1) ) 5 , 6 , 5 

5 CALL COUNT ( IT,N,3, 15,2) 

CONVERT Y TO SECONDS AND X TO MM FROM ORIGIN 
DO 204 1 = 1, N, 2 

201 XX= I T ( I ) 

YY = I T ( 1 + 1) 

11=11+1 

GO TO ( 202,203 ) , INDEX 

202 DX( I I ) = (XX-DXO)*SGN 

DY ( I I ) = Y SC AL UDYO-YY) 

ND= I I 
GO TO 204 

203 RX( I I ) = ( XX-DXO) *SGN 

CHECK FOR R BEFORE FIRST D AND ELIMINATE 
IF(RX(II)) 205,206,206 
205 1 1 = I I - I 

GO TO 204 

2 06 RY ( I I )=YSCAL*(DYO-YY)+DREF+DELAY 
5886 FORMAT ( • • ,D15.7) 

204 CONTINUE 
GO TO 4 

FOR THE FIRST LAYER, BOTH THE D-POINTS AND THE 
REFLECTION POINTS MUST BE READ IN. IF "INDEX" IS 
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C EQUAL TO ONE, THE D-POINTS HAVE NOT BEEN READ YET 

6 GO TO (61,65), INDEX 

ADD NUMBER OF SECOND PER SWEEP (SPS) TO 
DY IF DATA IS CONTINUED AT TOP OF RECORD 

61 DO 64 1=2,11 
CHECK = DY( I )-DY( I -1 ) +SPS/2.0 
IF (CHECK. GE. 0. 0) GO TO 64 

62 DO 63 J = I ,11 

63 DY ( J ) = DY ( J ) + SPS 

64 CONTINUE 
GO TO 3 

65 RETURN 
END 
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SUBROUTINE COUNT 



C COUNTS THE NUMBER OF NON-ZERO ENTRIES PER CARD 



SUBROUTINE COUNT ( K , M, I STAR , I STOP , I NT ) 
DIMENSION K ( 1 ) 

DO 2 I =1 STAR, I STOP , INT 

I F ( K ( I ) ) 2,1,2 

1 M= I— 1 
GO TO 3 

2 M=I STOP+INT-1 

3 RETURN 
END 
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SUBROUTINE INTERP 



THIS SUBROUTINE INTERPOLATES D'S FOR R'S. 

RX=U AND RY=S TRANSMITTED TO HONDO I BY COMMON 
ADD NUMBER OF SECONDS PER SWEEP (SPS) TO 
RY IF DATA IS CONTINUED AT TOP OF DATA 



SUBROUTINE INTERP ( P ERCE , A DENI , ADEN2 , I R EF , L , SP S » 
IMPLICIT REAL*8( A-H,0-Z) 

D I ME NS I ONTO ( 10) ,HH(10),V(10),W(10), SW(10) , 

ICW( 10) tTW(IO) ,D( 90), R( 90)»ZM< 90 ) , P ( 5 , 6 > , X < 5 ) 
DIMENSION RX ( 90),RY( 90 ) , DX < 2 5 ) , DY ( 25 ) 

COMMON TO,HH,V,W,TR, DR, SW,CW,TW,D, R, ZM , P , C , SVN , VH 
COMMON X,TIT,RX,RY, DX, DY , DX 0, D YO , Y SC AL , SGN 
COMMON M,IND,II,ND 

6 DO 9 1=2,11 

CHECK = RY(I)-RY( I-I ) +SPS/2.0 
IF (CHECK. GE. 0. 0 ) GO TO 9 

7 DO 8 J = I , I I 

8 RY( J ) = RY ( J ) + SPS 

9 CONTINUE 
60 J=1 

INTERPOLATE D FOR R. 

DO 65 1 = 1, I I 
IF(I-l) 50,52,50 

50 I F ( R X ( I ) -DX ( J)) 53,53,51 

51 I F ( ND -J) 53,53,52 

52 D X 1 = 0 X ( J) 

DY1=DY ( J) 

COMPUTE SLOPE OF D-LINE BETWEEN EACH POINT 
IN ORDER TO LINEARLY INTERPOLATE BETWEEN 
EACH POINT 
J=J+l 

DELXY= ( DY ( J ) -DY 1 ) / ( DX ( J ) -DX 1 ) 

DELXY IS SLOPE OF D-LINE BETWEEN J AND J+l POINTS. 
GO TO 50 

53 DD =DYH-(RX( I ) -DX 1 ) *DELXY 

RX ( I )=DD*PERCE 

RX ( I ) IS NOW THE SECONDS OF DIRECT TRAVE TIME 
CORRESPONDING TO R V ( I ) . 

IF NOT THE FIRST LAYER, CONVERT D'S TO KILOMETERS 
IF(L-l) 55,55,54 

54 RX ( I ) =RX ( I )*VH 

55 CONTINUE 

WRITE( 6, 101) ADEN1,ADEN2, I REF, II 
101 FORMAT ( 1 1 ' , ' STATION *,2A4,' LAYER NUMBER', 12, 
1T38 , ' NUMBER OF RE ADI NGS 1 , I 3/ / / ) 

RETURN 

END 
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SUBROUTINE LSFIT 



C SUBROUTINE TO FIT T/X DATA TO FOURTH-ORDER LEAST- 

C SQUARE POLYNOMIAL OR FITS T2/X2 DATA TO A LINEAR LEAST 

C -SQUARE LINE: 



SUBROUTINE LSFIT <LL,L) 

IMPLICIT REAL*8< A-H»0-Z) 

DIME NS I ONTO! 1 0 ) , HH ( 1 0 ) , V ( 1 0 ) , W ( 1 0 ) t SW( 10) , 
1CWI10) tTW(lO) ,D( 90) jR( 90),ZM( 90 ) , P ( 5 t 6 ) t X ( 5 ) 

DIMENSION RX ( 90),RY< 90 ) , DX ( 2 5 ) , D Y ( 25 ) 

COMMON TO,HH,V, W,TR,DR,SW,CW,TW,D, R, ZM,P,C , SVN,VH 
COMMON X,TIT,RX,RY,DX,DY,DXO,DYO,YSCAL,SGN 
COMMON M , I ND , KO 

C R CORRESPONDS TO REFLECTED TIMES 

C D CORRESPONDS TO REFLECTED TIMES 

IN=L L 
I N0= I N+ 1 
NQ= I N0+ 1 
KN= I NO- 1 
DO 1000 M= 1 , I NO 
DO 100 1. N=1,M 
P ( M , N ) =0 . 

DO 1001 1=1, KO 
I F ( D ( I ) ) 1001,9001,9001 

9001 IF(M+N-2)6000, 6000, 6001 

6000 TX= 1 . 

GO TO 6002 

6001 TX=D ( I ) ** ( M+N-2 ) 

6002 P ( M , N ) = P ( M , N ) + T X 

1001 CONTINUE 
P(M,NQ) =0. 

DO 1000 1=1, KO 
I F ( D ( I ) ) 1 000,9002,9002 

9002 I F ( M- 1 16003, 60 03, 6 004 

6003 TX=R ( I ) 

GO TO 6005 

6004 TX=D( I )**(M-1)*R(I ) 

6005 P( M,NQ) =P( M,NQ) +TX 
1000 CONTINUE 

DO 1010 N = 1 , I NO 
DO 1010 M= 1 , N 
1010 P ( M t N ) = P ( N , M ) 

DO 1022 K= 1 , KN 
DO 1002 1=1, K 
11=1+1 
KK = K+ 1 

A1=DABS( P( 1,1)) 

A2=DABS ( P ( KK» I ) ) 

IF ( A1-A2) 1020,1 021,1021 

1020 DO 1200 I I N= I , NQ 
AA= P ( K K , I IN) 

P( KK, I IN)=P( I , I I N) 

1200 P ( I , I I N ) = A A 

1021 IF(A1-1.D-30)1002, 1002, 1210 
1210 DO 1201 NM=I I , NQ 

1201 P(KK,NM) = P(KK,NM)-(P( I ,NM)*P(KK, I ) / P { I ,1 ) ) 

1002 CONTINUE 

1022 CONTINUE 

C SOLVE MATRIX 

DO 1005 K = 1 , I NO 
I=NQ-K 
11=1+1 

X( I ) = P ( I , NQ ) 

I F ( I - I NO ) 1006,1050,1050 
1006 DO 1060 K X= I I , INO 
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1060 X! I )=X( I )-(P( I ,KX)* X ( KX ) ) 

1050 A1 = DAB S ( P ( I , I ) ) 

IF (A 1-1. D-30) 150 0, 1500 1 1005 
1005 X< I >=X( I ) /P( 1,1) 

A= I N+ 1 

I F ( IN-1 ) 1701 , 1701,4 
1701 I F ( I ND ) 4 ,4 ,21 

21 DD0= 0 . 

0K=K0 

IF (L.FQ.1J GO TO 69 
WRITE! 6,696) 

696 FORMAT (' O' ,3X, ' ALL VALUES ARE REDUCED VALUES') 

WRITE! 6,501) 

501 FORMAT!' ' , 1 X , ' X2 • , 8X , • T2 • , 8 X , ' T2 ' , 1 IX , ' D I FF 

1 '/' ' ,T13, 'MEASURED' ,2X, 'COMPUTED 

1 '/' + ' , IX,8( '_' ) ,2! 2X,8( '_' ) ) , IX, 14! '_' ) /) 

GO TO 502 
69 WR I T E ( 6 , 500 ) 

50 0 FORMAT! 'O' , 1 X , ' X 2 ' ,8X,'T2* ,8X,'T2',11X,'DIFF 
1'/' ' ,T13, 'MEASURED' , 2X, 'COMPUTED 

!•/' + ', IX, 8(''),2(2X,8('_')), IX, 14 ('_')/) 

502 DO 1061 K=1,R0 
RR=0 . 

I F ( D ( K ) )9003, 9004, 9004 

9003 OK=OK- l . 

GO TO 1061 

9004 DO 1062 1=1, INO 

IF! 1-1 >6006,6006,6007 

6006 TX=X ( I ) 

GO TO 1062 

6007 TX=X(I )*(D!K)**I 1-1) ) 

1062 RR=RR+TX 

DK = R ( K ) -RR 
DD0=DD0+DK**2 

WRITE! 6, 1700) D( K ) , R ( K ) , RR , DK 
1700 FORMAT (3F10.6,D15.7) 

1061 CONTINUE 
ST U=OK 

COMPUTE STANDARD DEVIATION: 

STANDARD DEVIATION IS EQUAL TO THE SQUARE ROOT OF THE 
VARIANCE. 

DDO=DSQRT{ DDO/ ( STU-A) ) 

WRITE(6,1802)DD0 

1802 FORMAT! 'O' , ' STANDARD DE VI AT I ON ' , T 26 , ' = ' , D 1 5 . 7 
1 , T45 , ' SEC**2' ) 

STD=DDO 
VA= 1 . / X ( 2 ) 

WR I T E ( 6 , 1900 ) X ( 2 ) , X ( 1 ) 

1900 FORMAT!' SLOPE OF T2/X2 L I NF ' , T 26 , 1 = ' , F9 . 5 , ' 

1'/' ',' T-ZERO INTERCEPT* ,T26, ' = ',F9. 5) 

C SET STD EQUAL TO VARIANCE. 

STD=STD**2 
D1 1=0 . 

D2 1 = 0 . 

DO 6 M=1,K0 
I F ( D ( M ) )6,9010,9010 
9010 D2=D(M)**2 
D2 1 = D2 + D2 1 
D1 1 = D( MH-Dli 
6 CONTINUE 

D211=0K^D21-D11**2 
SX=STD/D2 1 1 

C COMPUTE THE VARIANCE OF THE INTERVAL SOUND SPEED. 

SVN=OK*SX 
SVN=SVN*VA**4 
SVN=DSQRT ( SVN) 

GO TO 4 

1500 W R I T E ( 6 , 1501)1 , I 

1501 FORMAT (2HP! , 13, 1H, , 13, 24H) =0, PROBLEM UNDETERMINED) 

4 RETURN 

END 
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SUBROUTINE TUPPER 



THIS SUBROUTINE IS USED TO REDUCE THE EFFECTS OF THE 
UPPER LAYERS BY COMPUTING THE TRAVEL TIME IN THESE 
LAYERS. THE COMPUTED TIMES ARE TRANSMITTED TO HONDO I 
WHERE THEY ARE SUBTRACTED FROM THE TOTAL TIMES TO 
OBTAIN THE REDUCED TIMES. 

NOTE THAT TI IS TRANSMITTED TO HONDO I AS TIT. 



SUBROUTINE TUPPER (L,ZZ,IK) 

IMPLICIT REAL*8( A-H,0-Z) 

DIME NS I ONTO! 10),HH(10),V(10),W<10), SW(10),X(10), 

1CW( 10) ,TW( 10) ,D( 90), R( 90),ZM( 90 ) , P< 5, 6 ) » XX ( 10) 
DIMENSION RX ( 90),RY( 90 ) , DX ( 2 5 ) , D Y ( 25 ) , Y ( 5 ) 

COMMON TO,HH, V,W,TR,OR, SW,CW,TW,D,R, ZM , P , C » S VN , VH 
COMMON Y , T I ,RX,RY,DX,DY,DXO,DYO,YSCAL,SGN 
COMMON M,IND,KO 
CHECK IF FIRST POINT 
IF( IKJ9990, 9990, 9991 

FIND THICKNESSES , COMPUTE SIN COS TAN 

9990 DO 1000 1=1, L 
SW( ] ) =DS I N ( W ( I )) 

CW( 1 ) = DCOS ( W ( I ) ) 

1000 TW ( I )=SW( I )/CW( I ) 

C=l. 

IF ( L-I ) 1003,1003,1002 

1002 L 1 = L- 1 

ww=o . 

DO 1001 1=1, LI 

1001 WW=WW+W(I) 

C=DCOS(WW) 

FIND ANGLE OF REFRACTION IN UPPER LAYER CORRESPONDING 
TO RAY REFLECTED NORMALLY FROM LOWER INTERFACE OF 
LAYER CONSIDERED 

1003 I = L 
Q=W( I ) 

IF (1.-1)11,11, 10 

10 Q=DARSIN( (V( I — 1 ) / V ( I ) ) *DS I N ( Q ) ) + W(I-l) 

1 = 1-1 

I F ( I - 1 ) 11,11,10 

FIND ANGLE OF INCIDENCE Q AT UPPER INTERFACE OF 1-1 
LAYER. 

11 Q =Q-W( I ) 

U=0 . 

TT = 0 . 

12 IF(I-L) 13, 14,14 
COMPUTE HH( I ) 

13 HC = HH( I )-U*DSIN( W( I ) ) 

SI X=DS I N( Q ) 

COX=DCOS( Q ) 

COMPUTE HORIZONTAL SEPARATION FOR 1-1 LAYER(U). 
U=U*CW(I) +HC * S I X /COX 

FIND ANGLE OF REFRACTION XW FOR NEXT LAYER UP FOR 
ANGLE OF INCIDENCE Q FROM LAYER BELOW 
TT=TT+HC/ ( V ( I i *COX ) 

XW=DARSIN< V< I+1)*SIX/V< I ) ) 

FIND ANGLE XW MAKES WITH UPPER INTERFACE Q. Q BECOME 
NEW ANGLE OF INCIDENCE ON UPPER INTERFACE. 

q= xw-w ( n i) 
i = i + i 
GO TO 12 

14 HC= ( TO ( I ) / 2 . -TT ) *V ( I ) 

C GET THICKNESS OF LAYER OF LAYER BEING CONSIDERED 

HH ( I )=HC+U*SW( I ) 

9991 1=0 

C COMPUTE ANGLE EACH RAY LEAVES SURFACE WITH RELATIVE TO 
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C VERTICAL 

XX ( 1 ) = ZZ+ W ( 1 ) 

924 1=1+1 

IF(L-I) 9252,9252,9250 
9250 SIX = DSIN(XX( I ) ) 

I F ( SIX-V( I)/V( I + l) ) 9253, 30, 30 
C COMPUTE THE ANGLE OF ENTRY INTO THE LAYER FOR EACH 

C POINT: 

9253 SIXW=V( I+1)*SIX/V( I) 

XW=DARSIN( SIXW ) 

XX( I + l )=XW+W( I + l ) 

GO TO 924 
9252 I = L+ 1 

X ( L » =X X C L ) 

926 1=1-1 

I F ( 1-1)21,21 , 928 

92 8 X( I-U =DARSIN(DSIN(X( I ) +W ( I ) ) * V ( I- 1 ) /V ( I) ) 

GO TO 926 

C MAIN COMPUTATION 

21 TT = 0 . 0 
U=0. 0 
1 = 0 

24 1=1+1 

HC = HH ( I )-U*SW( I ) 

S I X= DS I N ( X ( I ) ) 

CO X=DC OS ( X ( I ) ) 

HK--HC*( SIX/COX ) 

U=U*CW<I) +HK 

TK=HC / ( V ( I )*COX) 

TT=TT+TK 

I F ( L- I ) 252,252,24 
252 I =L+ 1 

HK=HK/ DCOS ( W ( L ) ) 

C COMPUTE REDUCTION TO REFLECTIONS TR AND SEPARATIONS DR 

TR=TK 
DR = HK 

T I =(HC/V(L))*2. 

PP=U 

26 1=1-1 

RR=HH ( I )-PP*TW( I ) 

CO X= DC OS ( X X ( I ) -W ( I ) ) 

HK=RR*DSIN( XX( I ) ) /COX 
PP = PP/C.W( I ) +HK 
TK=(RR*CW( I ) ) / ( COX*V ( I ) 1 
TT=TT+TK 

I F ( I-L) 262,261,261 

261 DR=DR+HK 
TR=TR+TK 

262 I F ( I - 1 ) 27,27,26 
27 TR=TT-TR 

DR=PP-DR*C 
3 RETURN 
30 DR=- 1. D20 
TR = 0 * 

T I =0 . 

GO TO 3 
END 
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SAMPLE DATA 



APPENDIX D. SAMPLE DATA 
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APPENDIX E. RESULTS 



STATION DIANNE NUMBER LAYERS 3 
TRIAL SPEEDS AND DIPS 



LAYER 


TRIAL TRIAL 

SPD. DIP 


1 

2 

3 


1.5000 0.0 

1.6000 0.0 
1.7000 0.0 


SONIC 


SOUNDING SPEED = 1.490 KM/SEC 
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STATION DIANNE 



LAYER NUMBER I 



NUMBER OF READINGS 46 



X2 



M Ia£U££Q £OMPlJI£D 



DIFF 



0.159342 
0.470273 
0.739195 
0.965480 
1.244668 
1 . 508562 
1.825360 
2.024092 
2.357128 
2.648469 
2.956782 
3.282065 
3. 546797 
3.781881 
4.065679 
4.359744 
4. 61 9971 
4.978674 
5.256501 
5.638670 
5.934096 
6.135233 
6.600064 
6. 919376 
7. 246231 
7. 524372 
7.807751 
8.038226 
8.44 96 23 
8. 689317 
9. 116848 
9.428513 
9.809425 
10.067556 
10.527353 
10.794704 
11.065407 
1 1 .408500 
1 1. 756831 
12. 181743 
12.325056 
12.614197 
12. 980338 
13.351717 
13.880447 
14.033398 

STANDARD 
SLOPE OF 



27. 

27. 

27. 

28. 
28. 
28. 
28. 
29. 
29 . 

29. 

30. 
30. 
30. 

30. 

31. 
31. 

31. 

32. 
32. 

32. 

33. 
33. 

33. 

34. 
34. 
34. 

34. 

35. 
35. 

35. 

36. 
36. 

36. 

37. 
37. 

37. 

38. 
38. 

38. 

39. 
39. 

39. 

40. 
40. 
40. 
41 . 



293261 
5991 18 
818630 
127409 
415657 
683037 
951670 
266656 
51 5341 
856153 
13021 1 
47 4544 
705186 
959897 
285595 
612999 
918542 
178226 
462718 
796208 
131402 
37 1869 
758426 
049804 
342435 
660864 
956103 
227839 
575206 
899287 
224837 
551857 
855026 
184872 
541734 
797678 
183226 
441344 
778199 
142608 
430125 
797577 
113893 
457981 
803539 
150566 



27.364615 
27.672436 
27.938668 
28. 162689 
28.439085 
28.700339 
29. 013968 
29.210712 
29. 540416 
29. 828843 
30. 134071 
30.456101 
30. 718184 
30.950917 
31.231876 
31. 523000 
31.780623 
3 2 . 13 57 3 8 
32.410785 
32.789131 
33.081603 
33. 280732 
33 . 740908 
34.057026 
34.380612 
34.655971 
34.936515 
35.164684 
35. 571966 
35.809262 
36.232516 
36. 541063 
36.918165 
37.173714 
37.628912 
37.893588 
38. 161583 
38.501244 
38 . 846091 
39.266752 
39.408632 
39.694881 
40.057359 
40.425022 
40.948463 
41.099885 



-0.7135450D-01 
-0.7331870D-01 
-0. 1200381D 00 
0. 35280460-01 
0. 2342762D-01 
0. 1 7301360-01 
6229822D-0 1 



>5944300-0 1 



- 0 , 

O.i 

-0. 2507517D-01 
0.27309710-01 
-0. 38606780-02 
0. 18442590-01 
-0. 1 299866D-0 1 
0. 8979450D-02 
0. 5371908D-01 
0 . 8999903D-0 1 
0.1379184D 00 
0. 4248 84 ID-01 
0. 5193265D-01 
0.7076253D-02 
0. 4979869D-0 1 
0.9113680D-01 
0. 1 .7517720-01 
-0. 7222040D-02 
-0 .38177070-01 
0.4893397D-02 
0. 1958790D-01 
0.631 542 8D — 0 1 
0. 32391690-02 
0.90024190-01 
-0.76792580-02 
0. 10793660-01 
-0.63139040-01 
0.1 1157250-01 
-0.87178410-01 
-0. 9590993 D- 01 
0.21642460-01 
-0.59900000-01 
-0.67892650-01 
-0. 1241445D 00 
0.214^2850-01 
0.10269670 00 
0.56534470-01 
0.329591 7D-0 1 
-0. 14492400 00 
0.50681820-01 



DEVI ATI ON 
T2/X2 LINE 



T-ZERO INTERCEPT 
INTERVAL SPEED 
T-ZERO REFLECTION TIME = 
PROPAGATED ERROR 
DIP OF LAYER 
DEPTH OF SEA FLOOR 



0.63899330-01 
0. 99000 
27.20687 
1.4900 
5.2160 

± 0.23287260-02 

0.0 

3885, 
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SEC**2 



KM/SEC 

SEC 

KM/SEC 

DEGREES 

METERS 
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STATION DIANNE 



LAYER NUMBER 2 



NUMBER OF READINGS 30 



X2 



VALUES 



ALL VALUES ARE REDUCED 

M£AiUR£D ZaMPUTfQ 



0.048589 
0.053972 
0. 058720 
0. 062745 
0.068639 
0.075737 
0. 084365 
0.091527 
0.101089 
0. 113960 
0. 127683 
0.141809 
0. 160847 
0. 180489 
0. 196513 
0.226165 
0.256824 
0. 29 2664 
0.334553 
0.391120 
0.415052 
0. 47 7440 
0.528239 
0.62 1924 
0.675188 
0. 748504 
0.883215 
1.000119 



0. 013895 
0. 008540 
0. 012488 
0. 020485 
0. 024050 
0.026166 
0. 021904 
0. 029197 
0. 035198 
0. 03 72 00 
0.050201 
0.050976 
0. 053300 
0. 059381 
0. 078500 
0. 084194 
0. 091382 
0. 104789 
0. 113446 
0. 125347 
0. 143980 
0. 161590 
0.181217 
0. 201968 
0. 221144 
0.248038 
0.262237 
0. 303986 



0.019739 
0.021427 
0.022915 
0. 0241 76 
0. 026024 
0.028249 
0.030953 
0.033198 
0.036195 
0.040230 
0.044531 
0.048959 
0.054925 
0.061084 
0.066106 
0.075401 
0.08501 1 
0.096245 
0.109376 
0.127107 
0.134608 
0. 154164 
0. 170087 
0. 199453 
0.216149 
0.239130 
0.281356 
0.318000 ■ 



STANDARD DEVIATION 
SLOPE OF T2/X2 LINE 
T-ZERO INTERCEPT 
INTERVAL SPEED 
T-ZERO REFLECTION TI 
PROPAGATED ERROR 
DIP OF LAYER 
THICKNESS OF LAYER 



-0.5843 72 9D-02 
-0. 1288677D-01 
-0.1 04 2 65 4 D- 01 
-0. 3691739 D- 02 
-0. 197345 8D-02 
-0. 2083175D-02 
-0. 9048658 D- 02 
— 0.4 001252D-02 
-0. 9960851 D-03 

-0.30302160-02 
0. 5669983D-02 
0.201655 6D-02 
-0. 1 62562 5D-02 
-0.1 7 029 2 OD— 02 
0.1 239395D-01 
0.87928250-02 
0.63707930-02" 
0.85436200-02 
0.4 070864D-02 
-0. 17595640-02 
0.9371501 D— 0 2 
0. 7425848D-02 
0.11 12972D-01 
0. 25145260-02 
0. 49953220-02 
0. 8908204D-02 
-0. 1911 887D-01 
-0. 14014210-01 



ME 

2 



= + 



0. 82838790-02 
0.31345 
0. 00451 
1 . 7861 
5.8053 

0.59694770-01 

0.0 

526.28 



SEC**2 

KM/SEC 

SEC 

KM/SEC 

DEGREES 

meters 
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STATION DIANNE LAYER NUMBER 3 NUMBER OF READINGS 32 



ALL VALUES ARE REDUCED VALUES 
X2 T2 T2 

M£MUE£D £I!E£L!IE£ 



OIFF 



0.027340 
0.025440 
0.026635 
0.030224 
0.035604 
0.045221 
0.055817 
0.070666 
0. 084917 
0.099064 
0.119355 
0.139051 
0. 154577 
0. 175825 
0. 194427 
0. 217424 
0.233800 
0.244261 
0.255025 
0. 265551 
0.268851 
0. 270837 
0. 266871 
0. 256989 
0.246093 
0. 231927 
0.214524 
0.189585 
0.155604 
0.123375 



0.002257 
0. 003936 
0. 003679 
0. 006926 
0.014095 
0. 010986 
0.015926 
0.022306 
0.027209 
0. 034296 
0. 036131 
0.041973 
0.052730 
0. 055410 
0.058062 
0.066329 
0.066879 
0. 079904 
0. 082052 
0.077297 
0. 089557 
0. 081146 
0.074851 
0. 072681 
0.072250 
0.072446 
0. 07 0957 
0. 054853 
0.039932 
0. 026457 



0.006382 
0.005787 
0.006161 
0.007284 
0.008968 
0.011978 
0.015294 
0. C19942 
0. C24402 
0.028830 
0.035180 
0.041344 
0.046204 
0.052854 
0.058675 
0.065873 
0.070998 
0.074272 
0.077641 
0.080935 
0.081968 
0. 082590 
0.081348 
0.078256 
0. 074845 
0.070412 
0.064965 
0.057160 
0.046525 
0.036438 



-0.4125055D-02 
-0.18513660-02 
— 0.2482090 D- 02 
-0. 35 80794D— 03 
0. 5126746D-02 
-0.9923697D-03 
0.63 10724D-03 
0.2364083D-02 
0. 28075970-02 
0.5 46627 2D-0 2 
0.9506 511 D-0 3 
0.62868050-03 
0.65267710-02 
0.25567200-02 
-0.61386410-03 
0.45612860-03 
-0.41194490-02 
0.56317190-02 
0.44111040-02 
-0.36382140-02 
0 .75890490-02 
-0.14438620-02 
-0.64975230-02 
-0. 55747730-02 
-0.25951320-02 
0.20346890-02 
0.59915750-02 
-0.23070990-02 
-0.65926600-02 
-0.9981321 D-0 2 



STANDARD DEVIATION 
SLOPE OF T2/X2 LINE 
T-ZERO INTERCEPT 
INTERVAL SPEED 
T-ZERO REFLECTION TIME = 
PROPAGATED ERROR 
OIP OF LAYER 

THICKNESS OF LAYER 3 = 



0.44653880-02 
0.31297 
-0.00217 
1.7875 
6. 2508 

± 0.952011 8D— 0 1 

0.0 

398. 16 



SEC**2 



KM/SEC 

SEC 

KM/SEC 

DEGREES 

METERS 
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SPEED OF SOUND AT WATER SURFACE = 



1.483 KM/SEC 
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